Linear Temperature Dependence of the Lower Critical Field Hd in F-Doped LaOFeAs 
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We present the first experimental results of the lower critical field Hd of the newly discovered F-doped 
superconductor LaOo.sFo.iFeAs (F-LaOFeAs) by global and local magnetization measurements. It is found 
that Hc\ showed an clear linear-r dependence down to a temperature of 2 K, indicative of an unconventional 
pairing symmetry with a nodal gap function. Based on the J-wave model, we estimated a maximum gap value 
Ao = 4.0 ± 0.6 meV, in consistent with the recent specific heat and point-contact tunneling measurements. 
Taking the demagnetization factor into account, the absolute value of //ci(0) is determined to be about 54 Oe, 
manifesting a low superfluid density for LaO0.9F0.1Fe As. 

PACS numbers: 74.20.Rp, 74.25.Ha, 74.70.Dd 



The recent discovery of novel superconductivity in rare- 
earth iron-based layered superconductors has received great 
attention in the scientific community III I2I [sl y, Isl la LZi 
[SIlElEllIillilElElQ- Except for the cuprate 



superconductors, these new materials ReOi_iF,iFeAs exhibit 
quite high critical temperatures: 26-28K for Re=La at x=0.05- 
O.llill M S H, 36K for Re=Gd at x=0.17|l(Jl, 41K for 
Re=Ce with x=0.16||l, 43K for Re=Sm with x=0.15|6], and 
52K for Re=Nd, Pr at x=0.15|li H, as well as 25K in hole 
doped Lai_vSrvOFeAs \5\. One of the crucial issues to un- 
derstand the underlying superconducting mechanism in such 
transition metal-based systems is the pairing interaction, i.e., 
the symmetry of the superconducting order parameter and 
the nature of the low energy excitations. Up to now, there 
have been several investigations on the pairing symmetry for 
LaOo.gFo.iFeAs (F-LaOFeAs) superconductor Recent spe- 
cific heat measurement has revealed a nonlinear magnetic field 
dependence of the electronic specific heat coefficient y{H) in 
the low temperature limit, which is consistent with the predic- 
tion for a nodal superconductor jlTT . The presence of a node 
in the superconducting gap has also been verified by the ob- 
servation of a zero-biased conductance peak in point contact 
tunneling spectroscopy 1 18]. On the other hand, an extended 
s-wave multi-band superconductivity was proposed jllll to ac- 
count for the experiment results of the extremely high upper 
critical field in F-LaOFeAs. From the band structure point of 
view, the two-band model has been strengthened with li-wave 
pairing interactions originated from the intra-band antiferro- 
magnetic coupling plus an effective inter-band antiferromag- 
netic interaction 11611 . Because the current experiments have 
been performed on polycrystalline samples, the experimental 
results within the context of the symmetry of the pair state is 
not yet well determined. 

Lower critical field Hc\{T), or magnetic penetration depth 
A(T) are fundamental probes of the existence of nodes, or two- 
gap in the superconducting gap function of unconventional su- 
perconductors. As an advantage, HdiJ) measurement probe 
relatively large distances {A ~ 1000 A) and are far less sen- 
sitive to sample surface quality. In this work we present the 
first detailed measurements of the temperature dependence of 
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FIG. 1: (Color online) Magnetization hysteresis loops of 
LaOo.cjFo.iFeAs measured by VSM in ZFC mode at T = 2, 10 K. 



the lower critical field Hc\ of superconducting LaOo.gFo.iFeAs 
polycrystals. We found a predominately linear-T dependence 
of Hc\ down to a temperature of 2 K, which is a strong ev- 
idence for nodes in the superconducting gap of F-LaOFeAs 
samples. 

Our polycrystalline samples were synthesized by using a 
two-step solid state reaction method. The superconducting 
transition temperature Tc defined as the onset of the drop in 
resistivity was 26 K with a transition width of AT^c = 3 K 
(10%-90% of normal state resistivity). The details of the sam- 
ple synthesis and the characterization have been reported in 
our previous papers isHlTtl- Our samples were relatively dense 
with a metallic shiny grain surface and shaped into various di- 
mensions. The data presented here were taken on a platelike 
polycrystal with the dimension of 0.6 x 0.6 x 0.2 mm^. 

Global dc magnetization measurement were carried out 
by a high resolution vibrating sample magnetometer (VSM) 
(Quantum Design) with applied field H parallel to the short- 
est lateral side of the sample. Local magnetization loops were 



measured using a two dimensional electron gas (2DEG) micro 
Hall probe with an active area of 10 x 10 jim^. The Hall probe 
was characterized without sample attachment at different tem- 
peratures. In our experiment, all M{H) curves were taken in 
zero field cool (ZFC) mode with initial temperature up to 40 
K. To minimize the complex effects of the character of the 
field penetration in layered structure, such as Bean-Livingston 
surface barriers and/or geometrical barriers 1. 19.. .20.] . we used 
a low field sweep rate of 30 Oe/min to measure isothermal 
magnetization curves. 

Typical global magnetization loops M(H) of our specimen 
at T-2 and 10 K are presented in Fig. 1. It is seen that 
the M(H) curves are symmetric and show a relatively large 
critical current density in the mixed state using Bean critical 
state model. This implies that the surface barriers for flux 
entry does not play an important role in our case since Bean- 
Livingston surface barriers are expected to give rise to a very 
asymmetric magnetization loop 12111 . and the pinning mech- 
anism mainly arises from bulk pinning. In addition, a fea- 
ture of a peak in M{H) appears followed the full penetration 
field. This so called "second peak" in M{H) loops has been 
widely observed in strongly layered superconductors, such as 
Bi2Sr2CaCu208 and NdCeCuO single crystals, and attributed 
to the vortex phase transition to vortex glass state 11221 12311 . 
Here the observation of second peak effect in our specimen 
indicates that the LaOo.gFo.iFeAs is strongly anisotropic and 
the intergrain links are strong, representative of the bulk prop- 
erty. 

Figure 2(a) and Figure 2(b) are the initial isothermal M(H) 
curves measured by VSM and Hall probe over the temperature 
range from 2 to 20 K, respectively. The second peak effect is 
more distinct in the curves by local measurement, as shown 
in Fig. 2(b). These M{H) curves show clearly a linear depen- 
dence of the magnetization on field caused by Meissner effect 
at low fields. In global measurement, due to the bulk pinning 
the entry of vortices is gradual, comparing to the local mea- 
surement. This causes a gradual M(H) curve, as shown in Fig. 
2(a). It is difficult to determine Hd using these M{H) curves. 
Therefore we use the method in Ref. 24 to extract the value of 
Hc\ : According to the Bean critical state model for type-II su- 
perconductors, magnetic induction B induced by the external 
field is: 



B = A(H - H.if/H* (H,i <H < //*), 



(1) 



where A is a sample-dependent constant and H* is related to 
critical current density. According to this model, a plot of 
B''^ vs H should yield a straight line with a threshold at Hd- 
Fig. 3 shows such plot of B^'^ vs // at T = 2, 3,4 and 10 K 
for the global curves, where B'^^ = (M - Mml)^'^ and Mml 
is the Meissner line for each curve, shown as the dotted line 
in Fig. 2 for the curve of T = 2 K. As displayed, the data 
are well described by Eq. (1) except a small kink at // > Hd . 
This kink is due to the second peak effect in the original M{H) 
curves. Neglecting this effect, Hd is well determined as the 
field of the threshold of B, as demonstrated in Fig. 3. By the 
same procedure, we also determined the Hd{T) data from the 
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FIG. 2: (Color online) The magnetization curve M{H) by VSM (a) 
and Hall probe (b), at various temperatures, respectively. The dotted 
red lines are the "Meissner line" showing the linearity of these curves 
at low fields. 



M(H) curves using Hall probe (not shown here for simplicity) 

Depicted in Fig. 4(a) are our main results, in which the ex- 
tracted Hd is plotted as a function of T for both VSM and 
Hall probe measurements. As shown in the inset of Fig. 4(a), 
the Hd measured by VSM is much lower than those by lo- 
cal measurement even though we measured the same sample. 
The reason for this difference is as follows. In global measure- 
ment, the VSM coil picks up any signal once vortices start to 
enter into the specimen through the edge of the sample. On 
the other hand, the local Hall probe detects signal only when 
vortices almost fully penetrates into the sample and reaches 
the region where the Hall probe locates. Assuming the Hall 
probe is located at the center of the sample, which is far from 
the sample edge, the Hall probe detects signal only from the 
central region. Accordingly, the Hd measured by local Hall 
sensor is expected to be larger than that of global VSM mea- 
surement. It also explains the experimental results that the 
M{H) curves by global measurement is more gradual around 
the penetration region, comparing those measured by local 
magnetization measurement, which is exactly the case we ob- 
served in Fig. 2. The fact that the demagnetization factor is 
different for local and global measurements can manifest itself 
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FIG. 3: (Color online) The M(H) data in the form of {M-MmlY'^ vs 
H. Mml is the Meissner line shown in Fig. 2. The curves (7=3,4,10 
K) are shifted up for clarity. The colored lines are guides to the eyes. 



by scaling the Hc\{T) curve measured by VSM to the one by 
Hall sensor. After we correct a scaling factor of 3.22 for the 
Hc\(T) by VSM, the two Hc\(T) lines almost collapse onto 
the same line except the part of high T, which is due to large 
error bars, as shown in the main panel of Fig. 4(a). The co- 
incidence of the two scaled Hc\{T) curves indicates that our 
results are intrinsic in physics, independent of measuring de- 
vice and technique. Thereafter, if not specially mentioned, we 
choose the data obtained by VSM for our analysis and discus- 
sions since the demagnetization factor is well determined for 
global measurement, as discussed below. 

We now analyze and discuss the obtained data of HdiJ). 
Principally one can get the absolute value of //ci(0) by ex- 
trapolating the nominal HdiJ) curve to T = 0, taking the 
demagnetization factor into account. For samples with ellip- 
tical cross sections, H^y can be deduced from the first pene- 
tration field //*j, assuming that the magnetization M - -H^ 
when the first vortex enters into the sample. Thus H has been 
rescaled to H - Ha - NM and Hd - H*^ /(I - A^), where A^ is 
the demagnetization factor and //„ the external field. For sam- 
ple with a rectangle cross section, the M{H) curves can still 
be approximately by a linear variation in the Meissner state. 
However, the demagnetization factor has been reconsidered. It 
has, for instance, been shown by Brandt that a bar with a rect- 
angle cross section, the first penetration field i/*j and Hd have 

the relation 1,25,1: Hc\ - //*j / tanh ( V036F7flj, where a and 
b are the width and the thickness of the samples, respectively. 
Using this formula, we estimate the effective demagnetization 
factor Neff - 0.67 as we take a ^ 0.6 mm and b - 0.2 mm 
for our sample. This estimated Neff is consistent with those 
of MgB2 platelets with the similar size ratio 126112711 . in which 
Neff - 0.6 - 0.69 has been determined. //*j(0) - 18 Oe is ob- 
tained by linear fitting to the data of Hc\{T) below 5 K. With 
Neff, the formula HdiQ) - //*j(0)/(l - Neff) yields a true 
Hc\(0) = 54 Oe. A striking feature of the nominal Hc\{T) 



shows an approximate linear-T dependence from around 6 K 
down to 2 K, the lowest temperature we can reach in this ex- 
periment, as shown by the straight line in the inset of Fig. 4(a). 
This approximate linear- T behavior is also evidenced from the 
Hc\{T) curve using Hall sensor. 

In order to investigate the properties of superconducting 
gap using the obtained linear-T dependence of Hd in the low- 
T region, we use the formula to calculate the absolute value 
of the penetration depth A (taking the demagnetization effect 
into account): 



<D, 



^--z7i2(i"'^+o4 



(2) 



where <I)o - hejle - 2.07 x 10"^ Oe cm^ is the flux quantum, 
and K is the Ginzberg-Laudau parameter Assuming k ~ 100, 
we deduce a value of A ^ 390 nm with the calculated param- 
eter //,i(0) = 54 Oe. 

In Fig. 4(b) we plot the deviation AA - A(T) - A(0) as 
a function of T using A(Q) - 390 nm. At low temperatures 
T < 5 K, AA shows a linear function of T, i. e. AA oc T, 
which is more explicit in the inset of Fig. 4(b). This lin- 
ear dependence of AA contradicts the result of an isotropic 
.s-wave gap superconductivity. For an isotropic .s-wave su- 
perconductor, such as conventional superconductor, the low 
excitation of the finite energy gap is a type of an exponen- 
tially thermal activation. In the isotropic .s-wave BCS theory, 
AA - A(0) ■\/7TAa/2T exp( - Ao/Tj, where Aq is the maximum 
gap value. This theory can not describe our AA{T) since the 
exponential term gives dA{T)/dT ^ at low T, shown as the 
fitting line in the inset of Fig. 4(b), in contrast to our result of 
dA(T)/dT oc T. 

The linear-T behavior of AA has been observed in high- 
Te superconductor YBCO Sill, BiiSraCaCuaOg jHISl 



single crystals, HgBa2Ca2Cu308+5 crystal powder 113 ill and 
magnetic superconductor YNi2B2C [32]. In cuprate super- 
conductors, such a linear-T dependence of penetration depth 
A at low r is a signature for the existence of a line node in the 
superconducting energy gap function, i.e. the li-wave pairing 
128: , . 3 1.1 . where thermal excitation of quasiparticles near nodes 
in the superconducting energy gap cause p, = 1 /A^ to decline 
linearly with temperature. Based on the li-wave pairing the- 
ory, Hci(T)/Hd(0), i.e., A^(0)/A^(T\ is a linear-r function at 
T «: Tc, and can be expressed as 1311 1330 : 



1 - 



HcijT) 
HdiO) 



1 



AHT) 



21n2— . (3) 

A() 



As shown in Fig. 4(a), the slope of the fitting line 
d(Hd(T)/Hei(0))/dT is 0.78, which yields a Ao = (1.8 + 
0.3)kBTc based on Eq. (3). Using the known Tc = 26 K, 
we found Ao = 4.0 + 0.6 meV. This deduced Ao are in excel- 
lent agreement with the values obtained by our recent specific 
heat 1 17] and point-contact tunneling spectroscopy [18] mea- 
surements on similar samples. 

Before conclusion, we would like to address the concern of 
the effect of impurity scattering on the linear-r dependence 
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FIG. 4: (Color online) (a) T dependence of H^i determined from the 
points deviating from the linearity of the initial M(H) curves by VSM 
and Hall Probe measurements. The data of Hci(T) by VSM is scaled 
by multiplying a number of 3.22 to account for the different demag- 
netization effects. The solid blue line is a rf-wave weak-coupling 
BCS fit with Ao = 3.85 meV; the dashed red line is an i-wave weak- 
coupling fit. Inset: the original data ofH^iiT) measured by VSM and 
Hall Probe, respectively. The solid lines are the linear fit to the data 
of low T. (b) r dependence of AA calculated with the data of H^i (T) 
and Eq. (2). Inset: the enlarged low temperature part, which shows 
a clear linear-J dependence of AA. The dotted blue line is the linear 
fit, and the dashed red line is an i-wave exponential-T line. 



of magnetic penetration depth. For our samples, the coher- 
ence length ^0 ~ 80 A 1,17.1 and the mean free path / > 100 
A yl], and thus our samples are in the moderate clean limit 
(/ ~ ^o)- Meanwhile, we believe that nonmagnetic impurity 
scattering does not strongly modify the low-T properties of 
F-LaOFeAs since Tc of our sample, another important low-T 
parameter, does not change much while impurity level (resid- 
ual resistivity) increases up to 2-4 times |3]. This situation 
also manifests itself by the fact that although different groups 
adopt different methods to synthesize F-LaOFeAs supercon- 
ductor with different impurity level, the Tc defined as the drop 
of resistance is around 25-28 K H fl H [11 [H . 

To summarize, we conduct both global and local 
magnetization measurements on the new superconductor 
LaOo.gFo.iFeAs. The temperature dependence of the lower 
critical field H^i is extracted. It is found that Hc\ show a 
linear-T behavior at low temperatures, suggesting a nodal 



gap function. We also obtained the lower critical field of 
Hci(Q) - 54 Oe, showing a low superfluid density for F- 
LaOFeAs. Based on the finite slope of Hc\{T) curve at low T, 
we estimated the maximum gap value to be A(0) - 4.0 + 0.6 
meV, which is very consistent with the results of recent spe- 
cific heat and tunneling experiments on the similar samples. 
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